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THE SEOUL AK PEETDEBATIONS OF TWO PLANETS MOVING IN 

THE SAME PLANE; WITH APPLICATION TO 

JUPITEE AND SATURN. 

By Db. G. W. Hill, Washington, D. C. 

The solution of this problem, when we restrict ourselves to the first 
powers of the eccentricities, is as old as Lagrange, and is well known. Le- 
verrier, in going over this ground, attempted to include the effect of the terms 
of three dimensions with respect to eccentricities and inclinations.* But when 
his method was applied to the four interior planets of the solar system it led 
to results that were nugatory. This method being that of successive approxi- 
mations, the expressions for the unknowns obtained in the simplest form of 
the investigation were substituted in the terms of three dimensions ; in conse- 
quence, he arrived at the same linear differential equations as before, but now 
augmented by known terms. His difficulty, in the case of the four interior 
planets, arose from the appearance in the results of integrating divisors which 
might receive very small, or even zero, values within the range of uncertainty 
of the values of the planetary masses. 

As far as the general question is concerned, no one has attempted to 
push the investigation further. Under these circumstances I have thought it 
might be well to treat as completely as we can the very simple case where we 
have only two planets executing their motions in the same plane. Although 
we see hero at a glance that the problem is reducible to quadratures, yet this 
taken by itself does not constitute a practical solution. Some difficulties are 
encountered in deriving from the quadratures series suitable for calculating 
the values of the unknowns. These difficulties I have succeeded in surmount- 
ing by a process which would not suggest itself, I think, at first sight. 

In the application which I have made to the case of Jupiter and Saturn 
with neglected mutual inclination, I have carried the approximation to quan- 

*AnnaIes de I'Dbservatoire de Paris, Tom. II, pp. 105-170 and pp. [38]-[51]. 
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tities of the fifth order, inclusive ; and it is not difficult to see what must be 
done if it is desired to go further. 

I. 

The first thing to be done in this investigation is to find a proper develop- 
ment of the potential or perturbative function. Quantities belonging to the 
interior planet will be denoted by symbols without an accent, and those 
belonging to the exterior by symbols having an accent. Let, then, m, r, a, g, 
u, and f denote severally the mass of the planet, the radius, the semi-axis 
major, the mean, eccentric, and true anomalies, while we denote the distance 
between the planets by A • The potential function Q is then given by the 
double definite integral 






or, if the integration is accomplished with reference to the eccentric anomalies, 
by the double definite integral 



in'' J oJo a a A 



These formuliB show that the potential function is proportional to the average 
value of the reciprocal of the distance when the mean anomalies are regarded 
as the independent variables, or to the average value of the product of the 
radii divided by the distance when the eccentric anomalies are the indepen- 
dent variables. As the eccentricities e and e' and the longitudes of the peri- 
helia 0) and w are the variable quantities whose forms as functions of the 
time we are seeking, it is plain they must be left indeterminate in the expres- 
sion we obtain for Q. Since A can be expressed in terms of u and u' as a 
finite form, the second formula for S2 is to be preferred. 

If Y be put for w — co', the expression for A > in the case we treat, is 

A = '-'[l - 2 5 cos (/-/' + r) + Ji]* • 

Thus, the expression for SJ becomes 

n 1 r^irZ-iTr r mm' , , , 

■J'S = -: — 2 I I , -z: sz^rr (lU au . 



''[l-2^,cos(/-/' + r)+^]^ 
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If Bj denote the same function of -, that Laplace's Sj'-'' is of a, the ratio of the 
mean distances, we may write 

[l - 2 J cos (/-/' H- r) + ^^ = \ 'j^^ B, cosj if-f + r) 
s denoting the base of natural logarithms. If we make s " «'— i = s, and put 



1 + t/I 



from the equations 



1 + V^ —e^ 



r = a (1 — e cos u) , r cos/'= a (cos u — e), r sin,/^ ai/1 — e^ sin u^ 

it is easy to derive 

r = a^^ (1 — cos) 



Thus 



s/v-i = 






5 — 10 
1 (U« ' 



S (0 


i 


s' — (y' ~ 


1 — (OS 


[l — co's' 



—J 



ewV— 1 . 



Seeking now an expression for Bj in terms of s and s', we have 

(1 — 2a cos 9? + a^)~i ■■ 



1 j = +» 



(we omit Laplace's subscript ^, as it is unnecessary for the purposes of dis- 
tinction). "We can regard 6"' as an approximate value of Bj , and the true 
value can be developed in a convergent series by Maclaurin's Theorem, if the 
perihelion radius of the exterior planet always exceeds the aphelion radius of 
the interior ; that is, if 



ae -\- ae 
a' — a 



<1. 
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The augmentation which « receives is 



rj{l — ios) 



-. — a = a 
r 



Thus 



"J — - ^T " -j-r 



Expanding the latter factor by the binomial theorem, 




v-1 






■q (1 — los) 



1-L" 






Substituting this value of Bj in the expression given ^.bove for — , and multi- 
plying the result by 

vin'ir mm' ,-, , 

r- = 1- '1 (1 — <"*') 

aa a ' ^ 

and employing the symbol V to denote the operation of taking the coefficient 
of sV in the development of a function of s and s' in a series of integral 
powers and products of s and s', we shall have 

^ mm'^=+» *=+"*=* (—1)'-* ,dW^ ,^. , ^ ■ .- 



jg-y) = 7/ V p (1 _ cosf-i \\ - jj '"^•'1 . 



X 



Let us put 



This quantity is then a function of e. Let E'^^'' be the same function of e' that 
^/* is of e. Then we can write 

mm' ■* =+"<= + " *=* i IV--* c/W-^' ■ — 
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This constitutes the infinite series to be employed in this investigation, 
and it remains only to study the properties of the functions of e denoted by 
7i7'"- l^y expanding the binomial factors involved in J?/-*' and performing the 
operation denoted by V, we shall get 

^«) ^ (_ ly (i+iH;l+2)^(i+i)./,„. 



XI 



. 1—.1 . 



(^'— ./ )^' — ■/ — 1) ^ { i — 1 ) 



1 i + 1 " 1.2 (j + 1) (J + 2) ' 

The series within the brackets is a case of the hypergeometric series 
1 4- «_ii^ + aja±l)lil-^ 1) a{a ^ 1) (« + 2)/3(/? + 1) (/9 + 2) , 



+ 



treated by Gauss in a memoir entitled " Disquisitiones generales circa seriem 
infinitam, etc." * This series gives the value of E^^'' in terms of tj and co , but 
it may readily be transformed into another expressed in terms of e. Adopt- 
ing Gauss's notation for this species of series 

^(.) = (_ 1). {i + ^){^+^^)...{ i+j) ,^^^>p^. _ . _ . . ^ i_ ^„.), 
But from Gauss's equation [100], p. 225 of the volume quoted, 



F{j - i, -i,J + 1, w') = (1 + wJ-^F 
and 



2 ' 



->i + i. 



(1 + < 



In consequence 



(1 + w' 



_ ( i + l)...(t+./) \_eyr-.^{i -./) {i ~J -1) 



ji -J) ji -,/ - 1 ) {i - J - 2) { i -J - 3) fej*. 
1 . 2 . C/ + 1) (.;• +-2) [2j + 



* See Gauss, Werke, Band III, p. 123. 
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It is remarkable that when i and _;' are integers the value of £'/-^' is equivalent 
to a rational function of the two quantities e and i/l — e^. For, when i is a 
positive integer, the series first given terminates after a finite number of terms. 
The same thing occurs in the second series when i — J is not negative. By 
Gauss's equation [82], p. 209 of the volume quoted, 



F 



'l±i .l+i + 1 / , 1 .2' 
2 ' 2 'J^'-'^ 



=(i-.^r^i^'ri=i±2,./=4+i,^+iy- 



From this it follows that 

_(^-l^..j^-i)^e 
~ 1.2....; 2 



{i -1) (i - 2) ... {i -j) U ^ ^ J _ ,2)-?^.f-> 



"j — M-2 j — i - 
' 2 ' 2 



",./ + l,e^ 



L' 



(1- 



„ 2i-ir 



1 I (i-i-l )( ^-i -2) 

({_—./ — 1) ...({ — / — 4) 
1.2".(i + l)0' + 2) 



which affords a finite expression for E^^' when i is negative. It will be noticed 
that E^^'> = 0, when i, not zero, is not greater than j. 

In order that the symmetry of the expression for Q may be seen, we will 
write the development of this quantity at length without the employment of 
the summatory signs ; 



la I 



(0) 



- a --^^- {ErE'r - ErE'^] 

-I- ^ a2 !L^ [EfE','^'» - 2E.f^E'^l + Ef>E'^l] 






2.3 do? 



+ 



SE.fE'^l + 3Ef^E'!^l — Ef^E'!^i] 

) 



»?.?/<-' (Same expression as above, except that 5, i?, and) 
^' j-fi" now take 1 as the upper index instead of O.j 

mm' (Same expression, except that b, E, and E'} _ „ 
«' (now take 2 as the upper index instead of O.j 

mm' (Same expression, except that b, E, and E'} „ 
a' (now take 3 as the upper index instead of O.j 
+ 
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It may be noticed that the terms in Ei''E'!^l,Ef'E'!^l, Ef'E'^, E^^E'^H, 
E^^^E'^l , EfE'!^l , etc. can be omitted in writing the expression, as the latter 
factors of these products vanish. However the symmetry is more apparent 
when they are retained. 

The following table exhibits the values of all the ^'s required in develop- 
ing iJ to the terms of the sixth order, inclusive. They are expressed as func- 
tions of e, and the finite form is given as perhaps more interesting than the 
development in ascending powers of e. 

^,«» = 1 ^0*"' = 1 

^ w = 1 + J e^ ^i»l = (1 - e')-i 

^-3(0) = 1 + 1^2 E% = {\-^)-i 

E^ = 1 + Se' + P' E'^ ^\}'^l ^'1 ^-^ ~ ^'^~^ 

Ef = l + 5/ + ^e* ^(?) = 1^1 + I e'l (1 - ^)-i 

Er = 1 + Y '' + f "' + il '' ^-' = [i + ^^ + i '*] (1 - o-^ 

Ef> = l + ^e^ + ~e* + ge« E<^1= [l + 56^ + ^«^](1 -eV^ 



E,^"=~e E, 



(1). 



1 - |/ 1 



J5y» = -^e ^i'>=0 



3 

2' 

■L 5 TJ^iW J- 



El'' = - I ^ - ^ «' E% = e{l~ ^)-i 

E-"' = — 3e — ? «' - ^ e' ^<" — T^ , _l ^ 






2 8 



^,- =. _ 4. - 15^ - ^1 e^ - j| .^ ^i^i = p . + 1^^ .' + j| /•] (1 - .0- V 



4 
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F,'^ = l{l-^/^- 


-l/l- 


-«T 


^e)_(l-l/l-^T 


2 


c 




^0 ■ ^2 


^/>=l«2 










E^I\=0 


F<-^> = 1^ 










E% = Q 


E^^ = ^e' + y* 










E% = \^{\-^ri 



Ef> = 7e'-+7e* + ^e^ 



E!?\=~eHl—e')-^ 



|2\-S 






(i-i/r=7y , (i-i/r^)^ 



e e 





_1(1- 


-l/l- 


-A 




4 


e' 




^3<^' = 


-5«» 






^4*'' = 


4 






^5*^' = 


— 14<3«- 


-V 




^/') = 


— 21e^- 


_p 








^„<^' 



^i'i = 



^n=o 



.£"1^ = 



.{i-i/r=?)' 






^i^^^L^'a-e^)-^ 



£n=p^+^*' 



5 5 

32 



■](l-e^)-v 



In the present investigation it will be more convenient to make use of a 
development of ^/-^ in powers of A \ ^ - \ = 0. By substituting in 



J 



the formula for E^^'' in terms of e the values 

2 



Jl-=0^(l-^«)f 
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making, for the sake of brevity, i — J = k, and carrying the development to 
terms of the sixth order, inclusive, we obtain 

^* -^ 1) ^ f ^ Lu/+i) ^i 



[ 



l.'2.(i + l)0- + -2) 



2"\.u+ir 2.4 J" 

r k(k — l){k — 'l){k--Z){k — 4:){k — 5) 
+ L 1.2.3.(i + l)(i-l 2)(i + 3) 
j + ^k{k — \){k — 1){k-^) 
~2 1.2.(JH l)(,; + 2) 

+ 2.4 l.(./ + l) 27476 ^J j 

Or, particularizing with respect to j, 

E^ = 1 -M' (^ - 1) ^^ + i (^• - 1) p''~y,^~^^ - 1] d' 

i ii _ 1) (i 2) (>: - 3) r ({ - 4) {I - o) -].« 

r(i-lHi-2K* -3)^:-4) 3 (i-l)(t-2) Li n .,) 
^L 1.2.2^3 2 O 274J )• 

^,. ^ ('•±.l)iM:.2) ^. |i _^ |-(izz2)(jzi3I_ 1] eA , 

»'<., _ (»' + 3) ( t + 2)(z + 3) .;, 
' ~ 1.2.3 

And, specializing still further, 

^,<»' = 1 ^/> = 1 

^7) --= 1 + 2d-' — W Ef.\ :-= 1 + 2PF + W -h 8/f 

^w = 1 + 6^2 _ g^4 ^(0,^ =, 1 ^1^ 6^2 + 2W + 80i?« 

^;») = 1 H^ 12/92 __ g^4 _ i2^« Et\ = 1 + 12(9^ + 78(9* A- 380^« 

J?,*"' = 1 + 1W + 10(9' — mff' E'l\ = 1 + 20^2 ._!_ 190(^4 _^ i260^'' 

^ w = 1 + 30/9^ -I 60fl< — 160^" E% = 1 -h 30(92 _y 390^1 j 33(5o,9« 

^^10) = 1 ^. 42(9^ + 168^' — 280(9« ^^'J, = 1 + 42^^ _|_ 714^4 _j , 7723^ 
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i:m = _ ^ r? + L|^ <y2 - ?|I ^*"j ^1') .^ tf r4 + 82^^ + ~~ t)''] 

E,'^" = - ^ [8 H 116(?^ H^ 596'] i:% = ^ fs -[ -^^ ff' + i^ l?*] 

^V'> = e-' [3 - ^^j i^.;,<2' = ^2 [1 + ^^j 

^5<"^ = <^^ [21 + 2i<y^] E^i\ = d^[^ + 39^^] 

^„») = ff^ [28 + 8419^] E^l\ ^ ff' [6 -f lOm'l 

E;''^ = d' [36 + 20id'] E'l\ = i!'^ [10 4-- 230(9^] 

E,^'^ = — 48' E„^'> = — t)' 

E,^'> = — 56d' E!!>^6' 



iiV" = — 120#-' E}l'l=--=10d-' 

Through multiplication we obtain 

^,«'>.£'V'" = 1 

E.f>E''^^ = 1 + 2i9^ + W — 2(?* + 4(?^^'^ + 4.6" + 0/?" — M'd" + Si?^*?'^ + 8^'" 
Es^'^ETl = 1 + 6(?2 + 6^ — 6^^ -f 36i9'(^'' + 24^'* 

+ ()i9" — 36i9^r + 144i9^^'^ -[ 80(?"' 
E^^)E'<^, = 1 + 12(?2 + 12(?'^ — 6i9* + 144^^19'' -f- 78^ 

— 12^" — 72(9*r -f 936^^2^ + 380r' 
^^(0)^"(o) ^ 1 _^ 20/?^ + 20r 4^ 10^' -f 400(^2^'^ + 190^'* 

— 60^'' + 200i9^r A- dSOOOW' 4 1260r 
E.^'^E''^, = 1 + 30#- + 30tf''' + GOd' + 900^^<y''' -f 390/^'-' 

— 1606" ^1800ti'6''- -+ mOOtl-d" + 3360^'" 
E^E"-^. ^ 1 -I 42^^ + 42r' + 168^y* H- 1764^-^'^ 

+- 714/^'^ — 280^/' + 7056^V/2 _^ 29988^^^" -^7728^'" 
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/^ HI 1 ff^lj'-i _J_ ^ H'^ 



4 '2 4 J 



. 10 — 25^^ — 95/'7'^ + ^^^ ^* — -~ ffW — -^^ «"' 
4 2 4 



167, 

ftW 



_ 4 _ 2f'y2 — 22^'^ + 1^ 0' — llff'lP — i|^ H 

] 

18 — 99(^2 - 261^'^ + -]^ 0* — -^ ^^r — ~~^- d'*] 
4 2 4 J 

■ 28 — 266/^2 — 574r + ~0* — 5453^^^ — i?5i^-r 

A 2 



] 



£;(i)^^'l» = dd' [— 40 — 580(?? — llOOr — 295(^* — 15950ff'd" — 15115^'*] 

^,® £" w = m" [3 — (?' + 3/^''] 
EPE'!^1 =: OW [15 — 5^2 + 135^] 
Eli)E"3, = ffW [63 H- 68i9^ + 819^'^] 
EPE'21 = 9W'^ [168 + 504(?2 + 2968^] 
Ej'^''>E'2l = ^'^" [360 + 2040(^2 + 8280^] 

^®^'1« = — 56^3(?'3 
^/»)^w = — 386i9-^^'3 
^/•^)^'W = — 1200(?»r 

If, in the expression for ii, we call the function of the eccentricities 

( i)t (py^ii 

which multiplies ^^ — r~ a* --j~r- in the coefficient of cos jj, J!//-", and A denot' 

ing the characteristic of finite differences with respect to the variable i, it will 
be seen that we have 

Then the expressions for Mp can be derived by considering the preceding 
expressions, taken alternately with the positive and negative sign, as the suc- 
cessive diflferences of these functions with respect to the index i ; and it will 
be advantageous to apply the process separately to each power and product 
of and 0'. The exhibition of this follows : — 
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Coefficients of cos 0;- ; 





Coefficients of d 





Coefficient 


s of ^2 a 


1 















0~ 


-1 


+ 1 






- 2 ~ 2 ^_ g 










0" 


-1-hl 


-1 

+ 1 


+2+f-6~ 


12 + 20 















« 0" 



n 


^-10 

V 2 



































Coefficients of dK 

























+ 2 + 


2- 6 












— 2 

— 6~ 

+ 6 + 


4 

*4 16 

12-18 


+ 16^^0- 60 


168 








0~ 


n+ 6 


+ 24 





















30 

12 



+ 42 



Coefficients of d'^O'-'. 



— 4 ~ ^ + 36 

+ 28 + ^2 _ 108 ~ 1^* + 400 Q^^ 
?« -U 9fifi — "00 

_ 48 - '"^ + 148 + ^^^ - 500 , «. , + 1764 

_L 79 ^ _ _ _ 9/14. . . . + 864 
+ 24 

24 



T ..V. ' — 500 

72 _ 96 — 244 3g^ 

Q-24 , o. +120 











Coefficients of b" 





4 

16 

18 
6 





"~ * + 24 
+ 20^54 



34 

24 

6 





58 

30 

6 



78 
+ 112 



36 



190 
200 
+ 124 



390 
324 



714 
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Coefficients of ^*'. 






*^ + 12 + •'■^ — 60 



160 



I 12 + 12 _ 36 ^ *^ + 100 ' loo — 280 

' 04. 4-52 — 120 

— 12 -^^ + 16 + ^^ — 20 

-20^ 8 + 48 + ^2 

+ 20 + 40 

Coefficients oi d'^d'^. 



+ 4+ 4— 36 

— 28"~ ^^+ 36+ '^2+ 200 .ac^(^ 

_l_ 4. I 272 — 1800 

— 24 ^^ * + 308 + ^'^ — 1600 , .op:c + 7056 
+ 288 + ^12 _,020- 1328 + 3656 + 52^6 

— 420 ~ '^^^ + 1308 + 2^2^ 
+ 180 + 600 

Coefficients of d'^d'*. 



— 8 "~ ^ + 144 

+ 128+ ^^6— 792 ~ ^^6+ 3800 ^.^^^ 
' 656 I 2864 — 11700 

— 528 ''*"' + 2072 + '''"'* — 7900 , , ^oqc + 29988 

I 14.16 5036 +18288 ' 

+ 888 + ^^^^ - 2964 ^^'^^ + 10388 

— 660 ~ 1^4^ + 2388 + ^^^2 
+ 180+ 840 

Coefficients of ^'''. 



— 8~ ^8+ 80 

+ 64 + ''^ - 300 ~ 380 ^260 oo«n 

— 164 "~ "^'^'^ + 580 "^ »«"_2100 , ,0^0 + 7728 
+ 188 + 352-640-1220^,268 + '''' 



- 100 - 288 + 408 + 1048 
-. -I- 19:0 



+ 20^120 
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Coefficients multiplying dd' cos ;- : 

Coefficients of ^°. 



40 



+ 2 


/\ 






+ 2 


0-4 






— 2 — 


4 + 6 + 10 


-18 




+ 


2_2- 8 


+ l0 + 2«_ 





0+ 2 


— 2" 


- VA 
























Coefficients of ^^ 


■ 


Z\ 


— 2 






- 3- 2 


+ 23 + 25. 


— 99 




+ 18 + 21 
-12-30 


-51-^4. 
+ 42 + 93. 


+ 167 
-147 


+ 266 
— 314 


+ 12 


-12-54 

















Coefficients of d'^ 


. 


+ 1 








+ 1 '- 


-22 






-21-22. 


+ 73+ 95. 


-261 




+ 30 ^ 51 - 
-12-42 


- 93 - 1^^ 

f54 + 14^ 


+ 313 
-213 


+ 574 
— 526 


+ 12- 


-12- 66 







580 



— 1100 



Coefficients of ^* multiplied by 4. 



— 1 




- 1 0_^ 26 




+ 25+ 26 _ 159- 185 _^ ,.^3 




— 108 — 133 + 169 + 328 _ gi 


- 574 


_ 72+ 36^ 436+ ^67 _ jg^g 


-1754 


+ 400 + 472 _ JJJ2 - 1548 




- 240 - 640 





1180 
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Coefficients of O^'d'^ multiplied by 2. 



31900 



- 1 0__ 22 








— 23^ ^^ + 453 ' 


475 


■ 2871 


^ /% /^ /% /^ 


+ 408 + "f^ — 1943 ~ 


-2396 


8035 ' 


1- 10906 


— 1104 ~ -^^^^ + 3696 '' 


- 5639 


12959 


- 20994 


+ 1080 + ^^^^ — 3624 " 


-7320 






- 360-1***^ 









+ 3 

_^ 3 ^—334 



Coefficients of d"^ multiplied by 4. 



- 331 ~ ^^^ + 1861 + ^^^^ - 8451 , „,.ao 

+ 1196 + 1^27 _ ,3,5 - 6256 ,^231 + If 'I - 60460 

_ ,6,2 - 2868 55go + '''' - 19547 " ''''' 

+ 1040 + 2712 _ 3992 - 9572 

_ 240-1280 

Coefficients multiplying d^O'^ cos 2^ : 

Coefficients of ^"'. 

+ 3 

+ 3^0 

+ 3 ^-15^'^^+ 63_,gg 
-12-lJ + 33 + ^8_iog 3go 

+ 6 + 1^-24-^^+ 87 + 
0- 6+ 6 + 30 











Coefficients of ^^. 



1 

1 

, ^ + 5 



^- 



1 "+ 5"^ " + 63 

4 + ^ + 73 + ^^ - 441 

82+ 78-300 -3^3 + 1095 + 1^^^ 
,,Q — 222 .-- ^ 792 

60 + 200 
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Coefficients of t)'^. 

+ 3 

+ 3 » " 



3 _ J3J5 - 135 g^g 

— 1 Sf) -I- fi84. — 2968 

132 ^'^•^+549+ ^'"*-2149 , ^oio+8280 

I 4.1 4. 1 4.fi5 + 5dl2 

282 + *^* — 916 ^*'^'' + 3163 

220-^02+782+1698 

Coefficients of <9'^'^ cos Sy. 



+ 60+280 



+ 4 

+ 4 

^ ,. .„ 
+ 4 — 56 , QQc 

+ 4 " — 56 ^^ + 280 Q. , •- 1200 
-52- 56 + 168 + 224 _, 34 -864 

+ 60 + 112 _ 192 - 360 

-20- 80 



We can now write the explicit development of iJ as follows ; 



mm 2 i 



« ,, J. T ^y„) 



+ a ~ m' 4 2(9'^ — 2<9* -I- 4^^(?'^ + 4^ 4 Ol?" — 4:0*6'^ -\- SdW + S/?'"] 

-f J a^^ [2(?2 4- 2<9'^ — 2^^ + ^SffW + 16(?'^ 

^ «« + Qd^ — 28d*d'' + 128^2^'* + 64<9"'J 

+ 1 asf^°_' [6^< + 486l^r^ 4 18r — 12d' 

1.6 da ^ 24(9^r 4 528^^^'^ + 164<9'«] 

+ 1^ «<^' [6(?< 4 24^^r + 6r — 12<9'' 

^•^•* ^" 4 288^*1?'^ 4 888^^1?'^ 4 188^] 

+ ^-^^ a-^^'"' [201?'' 4- 4206l^(?'^ + 660<92r 4 lOOl?'"] 
2.3.4.5 ^a'^ '- ■' 

4- 1 a ""^ [201?'' 4 180«^^'^ 4 180(9^^ + 20^'"]! 
2.0.4.5.6 aa" t 
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+ 



2 df/r [_ 



"4 ~^ 2 



331^,/ 



_ 1 tf'.'^' [18^2 ^ 30^-2 _ 27(/^ + 204cdW + 299(?"] 
2.0 ata' 

_ _4-r«'^?r [12/^' + 12(?'^ + 18(?* + 552^2(?'2 + 4:186"] 
2.3.4: da* ^ ' -' 

— s-A— r «' ^* [100^* H- B^Od'd" + 260^'^] 
2.3.4.5 <^6f' '- ■' 

~273:"4T5.'6""S^ t^^^' + ^^^^'^" + *^^^"^ } ^^' '°' '' 



+ 



^^a^ + l^^]^-e^ + se^] 



da 



da' 



J 



2.3 « <r 



1 .747,(2) 



' 2.3.4 ri«' 
+ 2-73^6""^' t60^^ + 60^'^]| eW' cos 2^ 



, (^7.8) ^ ^*® , o 2dW> 2 ,(^»i('" 13 .(Z^Jf'* 



ofa 



<^«2 3 <^r/ 6 da: 

1 ..e^"*®) 



— Za"^ "^" ^-. «« "^" (93/9'3 cos 3r 

2 rf«-^ 36 <?«")■ ^'^ ''*'^'^'^- 

In order to have as few functions of a to deal with as possible, we gather 
together all the terms having the same powers of d and 0' as factors. Also 
it will serve our purposes better to have the development of ii in powers 
of cos Y than in cosines of multiples of ;-. For convenience in writing we 

denote a' --^-j- by (_;', i). We then put 



Ar=i0,l) + ~{0,2), 

A^ = - (0,1) - 1 (0,2) + I (0,3) + I (0,4) , 
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A,^"^ = 2(0,1) -h 7(0,2) + 4(0,3) + | (0,4) - 3(2,0) + 3(2,1) 

-|(2,2)-2(2,3)— J(2,4), 
Af^ = 2(0,1) + 4(0,2) H- 1 (0,3) + ^ (0,4) , 
Ar = - (0,3) - 1 (0,4) + ^2 (0,5) + ^2 (0,6) , 
^/«) = - 2(0,1) - 7(0,2) -I- 2(0,3) + 6(0,4) + ^ (0,5) + J (0,6) 

+ (2,0) - (2,1) -[- 1 (2,2) + |(2,3) - g (2,4) - ^ (2,5) - ^ (2,6) , 
Ar = HO,l) 4 32(0,2) + 44(0,3) + ^(0,4) + ^^ (0,5) + ^(0,6) 

- 3(2,0) + 3(2,1) - 1 (2,2) - 22(2,3) - ij- (2,4) - ^ (2,5) - A (2,6) , 
^,<») ^ 4(0,1) + 16(0,2) + ^ (0,3) + i^ (0,4) + ^^ (0,5) + i- (0,6) , 

A«- 2(1,0) -2(1,1) -(1,2), 

^« = - (1,0) + (1,1) -| (1,2) -3(1,3) -|(1,4), 

J, <■) = (1,0)- (1,1) -|: (1,2) -5 (1,3)- 1 (1,4), 

^3C) = -l(l,0) + i(l,l) + f(l,2) + |(l,3)-|(l,4)-|(l,5)-^(l,6), 

^«=-|(l,0) + |(l,l)-|^(l,2)-34(l,3)-23(l,4)-|(l,5)-|(l,6) 

- 12(3,0) + 12(3,1) - 6(3,2) + 2(3,3) + ^^(3,4) + 1(3,5) + ^ (3,6) , 

^/« = |{l,0)-|(l,l)-^|^(l,2)-2|?(i,3)_^(l,4)_l|(l,5)_,^l_(l,6), 

^/« = 3(2,0) - 3(2,1) + 1 (2,2) + 2(2,3) + ^ (2,4) , 

A\^'' = - (2,0) + (2,1) - 1 (2,2) - 1 (2,3) + g (2,4) + ^ (2,5) + ^ (2,6) , 

A/-^ = 3(2,0) - 3(2,1) + 1 (2,2) + 22(2,3) + ^ (2,4) + ^ (2,5) + ^ (2,6) , 

A,^^> = 4(3,0) - 4(3,1) + 2(3,2) - -| (3,3) - 1| (3,4) - 1 (3,5) - 1 (3,6) . 
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Then, neglecting the term which is independent of 6, d', and ;' for the 
reason that it is useless for our purposes, we shall have 

— , SJ = A,<-'>\d' + d") + A.'-'W + A.fWd" + A^V* + Afd" 

+ A^'Wd" + A^'-Wd" + A^W 
+ [^„(« + A^W + A,^W +A,<W + A^<Wd" + A.^'W'] 66' cos y 
+ [^o"' + -^i*''^' + A^W] 6'd" cos^ r 
+ A,^'>6W' cos' r ■ 

In order to make an application of the method to the case of Jupi- 
ter and Saturn, we take from Runkle's Tables of the Coefficients of the 
Perturbative Function the values of log (_;, i) corresponding to the argument 
log a = 9.7367414. 



i. 


j = 0. 


i = i- 


i = 2. 


i = 3. 





0.38^5227 


9.7929622 


9.4112303 


9.0721143 


1 


9.6447549 


9.9080135 


9.7808244 


9.5982418 


2 


9.9328686 


9.8807530 


0.0203420 


0.0219693 


3 


0.2943862 


0.3204279 


0.3188228 


0.3995660 


4 


0.8787099 


0.8712079 


0.8884960 


0.9011936 


5 


1.5571487 


1.5610571 


1.5658243 


1.5798073 


6 


2.8402885 


2.3412199 


2.3462289 


2.3533961 



Making use of these values, we obtain for this special case 

— , a = 0.8692176(<92 + 6'^) + 1.05019^* + 11.85269^^^'^ + 8.14486<9'* 

mm ^ ' 

+ 2.207^« + i6.1266*6'^ + 157.4:64:6^6'* + 89.730(?'« 

— [1.1365062 + 10.94248<92 + 22.34085(?'2 + 42.355(?* 

+ 335.261(?2<9'2 + 362.933^'*] 66' cos y 
+ [6.63740 + 86.288(?2 + 223.228^'^] 6^6'^ cos^ 

— 43.209<93^'3 cosV- 

II. 

The portion of the subject which treats of the integration of certain dif- 
ferential equations is now to be attended to. Denoting the mass of the sun 
by M, and putting 



fJL = M- 



M- 



G ^^m. \/ fm i/l — 6^, G' = m! \/ fia! i/l 
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the diiferential equations which determine the eccentricities and positions of 
the perihelia of the two planets are 



dG 


dQ 


dm dQ 


dt ' 


dm' 


dt dG' 


dG' _ 


dQ 


dm' dQ 


dt 


~ dm' ' 


dt dG' ' 


and 


co' only 


through Y =^ CO — 




dii 
dm 


dm 



Hence 

G -\- G' = a, constant 

is an integral of the problem. This integral equation may be more suitably 
expressed in terms of the variables d and 0' which we have before employed. 
Then K denoting an arbitrary constant, and denoting the constant quantities 

/— , -r-, V. 1 1 
m i/fxa, m V /xa by p , p , 

p + r^- - ^- 

The value of K is ascertained by substituting in the left member of this 
equation for d and d' the values they have at a definite epoch. We can 
now reduce the number of variables in the problem from four to three by 
adopting a variable v to replace d and 0', such that 

d = A \/K sin ^v , 6' = /' \/^ cos ^v . 

^v remains always in the first quadrant. Denoting the angles of the eccen- 
tricities by (f and f', the eccentricities are determined by the formulae 

e = sin <f, e' = sin f, 

sin ^(f ^ X ]/J{ sin ^v, sin ^f' =r /' ylC cos ^v. 

Making the substitutions in Q necessary to make it involve v instead of d and 
§', we put 

0^ = iPK{l — cosi^), ^'2=i;/2A"(l ^cosv), dd' = iU'Ksmu. 
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The function ii becomes, then, divisible by -ff", and, in order to simplify, we 
shall put a = Kli. Therefore, if we write x for cos v and put 

CI i Ji 4 



+ \ (/"J f'> + X'PA^<^ + ^^^'M/o) + A''Af'>) K^\ , 






i?,«» = 'I!^^ |1 {A' A i'' — n"Af^ + /'M/«>) ^ 

A<" = ^' ^^ |yl„'" + \ {A' A /» + a" A /") A^+ -J (^M/^ + )?PA^'^ + A'M,») ^^1 , 
<X A i Ji 4 L 

y^^o) == ??i^' .^Jjl (/^ yi 3O) _ n"AP + A'M,,<") ^^l , 

a o 

» ~" a' 8 " ' 

we shall then have 

+ [5;'> + /A<"« + /^/"a'^' + . . .] sin V cos r 
-y. [y,> (2) + y? P)^ _!_...] sin^ V cos^" y 
+ [^0*'" -f- • • •] sin^ V cos-'* ;- 
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With this expression for R it is readily seen from the preceding differen- 
tial equations that the differential equation determining v is 



dv 




1 


dR 


dt 




sin V 


dr' 




dx 


dli 






di' 


- dr 


' 



or 



Since R = a, constant is evidently an integral of the problem, we shall have 

dR dv dR dy „ 

dv dt dy dt 

Whence is derived 

dy ^ 1 dR 

dt sin V dv 

We still need an additional equation giving the value of some other function 
of Si and lo' than w — m . If we select <3 -\- Co' we have 

d {Co + Co) dii dii 

m da~dG'' 

If K is kept evident in the expressions for the various -5's, so that the par- 
tial derivatives of them with respect to this quantity may be taken, we shall 
have 



dR 



dQ _ d{KR) dK j^dR^dv^_ \d(KR) 

dG dK dG^ dv dG~ 2 dK „ , v dv ' 

2 tan jj 

dii __ djKR) dK irdRdv ___! djKR) 1 , vdR 
dG' ~~ dK dG' '^ dv dG' ~~ 2 dK ^ 2 2 dv ' 

Whence 

d {(0 -f- (5') d (KR) cos v dR 

dt dK sin v dv 

^diKR) dy 

In making our numerical application we take the mean distance a' as the 
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unit, when a becomes the same as « previously given, and assume for the 
masses the values 

1 , 1 

"1047.879' 3482.2" 

These give 

log X = 1.5758667, log A' = 1.7708956. 

The values adopted for the eccentricities at the beginning of 1850 are 

e = 0.04825801, e' = 0.05606467. 

These furnish the equations 

= [8.4778154] sin ^ , d' = [8.6728444] cos ^ , 

and the function i? becomes 

H = 0.0005906465 + 0.0002543964* + 0.00000196780*^ 
+ 0.000000019894«3 

— [0.0003548741 + 0.00000629406« + 0.00000008731«2] sin v cos y 
+ [0.00000148778 + 0.00000004479«] sin^v cos^ 

— 0.000000006560 sin'v cos'V- 

The value of the constant in the integral equation 

Ji=C 

is ascertained by substituting in the expression for 7i the values which u and y 
have at a definite epoch, as 1850. C being determined, the equation 11^ O 
can be solved, regarding sin v cos y as the quantity whose value is to be 
obtained. This value can be supposed developed in powers of cos v ^x, 
and we write 

sin V cos y = H = D^ + DyX + -^2^^ + ^i^^ + . . . . 

The readiest method of obtaining the Z>'s is by substituting the last expres- 
sion in R and then equating the resulting coefficients of each power of x to 
zero. We thus have a system of equations determining the Z>'s. These can 
be solved by successive approximation. If C is allowed to appear as an inde- 
terminate in the expressions for the !>'&, JI can be partially differentiated with 
reference to this quantity. 
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We can now make H play the role of R ; for we have 

dx dR dr dR , II 

- — -^ — — and y = arc cos - 



dt~ dr' dt~ dx' ' ~ yr-x^ ' 

Thus 

dH 

dt dy dC 

dx~'dC~ i/i_a,.2 —H^' 

where the radical in the denominator must receive the sign of sin ;' ; for we 
have 

cos V = £» , 
sin V cos y = 11^ D^ -\- D^x + Z>2«^ + DgX^ + . . . , 
sin y sin 7- = i/(l — x^ — 11^) . 

If we suppose the orbits are always ellipses x cannot pass the limits ± 1. 
Thus X must oscillate between a maximum and a minimum value, while 
dll/dC remains constantly of the same sign. The maximum and minimum 
values of x are evidently the two consecutive real roots of the equation in x 

\-x^—H^ = 0, 

which contain between them at any time the actual value of x. Calling these 
roots a and b, we may write 

1 — «2 — ^2 = (a — x) {x — b) Q, 

where Q is positive for all values of x lying between a and b ; and when the 
eccentricities are always small, the variation of Q is slight in comparison with 
its magnitude. In the place of x we can adopt a new variable, (p, such that 

a + b a — b , 

Then 

dx J I 

, = dip, 

l/(a— «)(« — b) ^ 

and the differential equation giving (p in terms of t is 

dR 

d£ ^ dC 

di'^ VQ' 
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To see how all this applies in the case of Jupiter and Saturn we assume 
the following values of the longitudes of the perihelia at the epoch 1850.0 : 

oj — 11° 54' 31".18, <ry = 90° 6' 57".55. 

The value of the constant C being now determined, and the equation li =^ C 
modified in such a way that it becomes more suitable for solution, we have 

0.4021256 = — 0.7168638* — 0.0055451*2 - 0.0000546*3 
+ [1 + 0.0177360* f 0.0002460*8 1 sin V cos ;- 
— [0.0041924 + 0.0001262*] sin^v cos^j- 
+ 0.0000185 sin^v cos';-- 

When this equation is solved with reference to sin i' cos ;- as the unknown, we 
obtain 

//= 0.4028046 4- 0.7121389* — 0.0050141*2 _ 0.0000050*^ 

And when we ascertain what increment // receives from an infinitesimal incre- 
ment in the quantity C, it results that 

_ 4{{= 2827.425 — 33.179* + 0.005*^. 

The equation 1 — x^ — ^^ __ q^ in this case, is 

0.8377485 - 0.5737057* — 1.5031028*^ + 0.0071447*=* — 0.0000180** »= 0. 

The consecutive real roots of this which contain between them the value of * 
at 1850.0 are 

a = 0.5803236, b = — 0.9586738. 

We derive from these the limiting values of v, which are 

54° 31' 36".14 and 163° 28' 14".01. 

Thus, when ;' ^ 0°, the minimum e of Jupiter has place, which is 0,02752623 ; 
as also the maximum e of Saturn, which is 0.08362800. And, when j- — 180"', 
the maximum e of Jupiter has place, and is 0.05944555 ; and the minimum e 
of Saturn, which is 0.01353514. 
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The remaining factor of the equation, two of whose roots we have just 
obtained, is 

Q = 1.5058180 — 0.0071522,r ^ 0.0000180a;2. 

"Whence 

~ = 0.8149177 + 0.0019353a; + 0.0000020«8. 
vQ 

Substituting, then, for x the expression 

X =- — 0.1891751 — 0.7694987 cos ,/', 



we get 

=r- 2304.1185 — 21.5662* — 0.0543* 



dt „on. ..or o. r^^o_, nnr.o 2 



= 2308.1802 + 16.5794 cos </' — 0.0161 cos 2</'. 

Integrating this, e being the arbitrary constant, 

t + c = 2308.1802^'- -I 16.5794 sin </> — 0.0080 sin 2i/'. 

Inverting this series and changing the numerical coefficients into seconds of 
arc we get 

</' = 19".05825 (i + c) — 1481".57 sin [19".0582o {f + c)] 
+ 6".04 sin 2 [19".05825 {t + c)]. 

From the value which (/' must have at the epoch 1850.0, t being counted 
thence, 

19".05825c = 277° 9' 9".15. 
Also, we have 

cos V = — 0.1891751 — 0.7694987 cos </', 
sin v cos r =^ + 0.2679063 — 0.5494490 cos </> — 0.0029673 cos^^'- 
+ 0.0000023 cos^'-, 

sin V sin y = [0.9446898 4 0.0017265 cos </> + 0.0000018 cos^^/-] sin </'. 

These equations enable us to determine the eccentricities and difference of the 
longitudes of the perihelia at any given time. 



Or 



Or 
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It remains to find the longitudes of the perihelia themselves. We have 
d(b 1 /-f 1 1 jT-dli L -\- X dy 

dt-^^^~^i^dK^ ^^>r dt ' 

dm _L /7 I ' TT^dll 1 — xdy 

'di '~ ^ ^ ^1 cfK ^~ Tt • 



d [m — ^ Ct) 1 T^ dy 1 i- X dy 

dx "'^^ dK ^ a~ dx ' 

' dx' ^ 2 dK 2" dx ' 

d (o> - \ Ct) .__ 1 dK ~^^ ''''''lix~ 1^3 



dx "2 1 (1 -^x" — IP) 

j^ dll .^ < dll Hx 

dim' —J Ct) ^ 1 dK^ ^ ^' Tlx ^ r+ ~^ 
dx "' 2" , (1— «2 — 7/2) 

Here K must be left indeterminate in the coefficients 2),, , Z*, , etc. of H, in order 
that we may get -^^ . In the next place, we derive 

K'^-n+x)"^^^- ^*' 
d (w — 1 (^) ^ 1 ^'- ^^r ^i h ^; ^^^- ^ _ ^ 

^i/^ " 2 ' ^ ^ Q 

j.dll , , _..\dll 4 Hx 

d [o/ - 1 C t) ^1 dK ~ ^ ^ dx ^ 1+x 

d4> 2 1 (? 

AVhen H, which is an infinite series in integral powers of x, is divided by 1 — x 
or 1 + «, remainders independent of x are left over which are equivalent to 
what H becomes when in it we make x ^1 and x = — 1. These remainders 
we denote as //(I) and //( — 1). Then we may write 



d{(o — ^^Ct) 



d4' 2 1 (> 

d4< 2 1 Q 
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Tlie difference of these equations gives 

17/(1) l7/(-l) 'T.T :iy in .nl,/ 

Since ;- returns to the same vahie after 9'' has augmented by a circumferenci', 
it follows that when the right member is expanded in an infinite series con- 
taining, besides two terms in the form of fractions having 1 — x and 1 -]- x 
as denominators, a set of terms proceeding according to cosines of multiples 
of f, the coefficient of the zero multiple of {.'■ must vanish. This is not imme- 
diately evident from the form of the expression. Hence I proceed to prove 
it to the degree of approximation we adopt. Let 

I Q 

then, omitting the two terms in the form of fractions and having 1 — ,'■ and 
1 -\- X for denominators, it will be perceived that we have 

'^^ = I),E„ + (2>„ -h A) ^, H I^^K - [AA; - A,E,\ X - [2 A£;, -^ n,E,\xK 

Substituting for x its value in terms of (j', if our proposition is true we ought 
to have 

AA; + (Z>„ + A) ^>; -f D.E, - \D,E„ - DJQ ""^^1 

-3 fa + b 1 - 1 

But if 



- [2 A A'„ - I),E,] g K+-^] - \ ab] = 0. 
Q = M„^ M,x H M.pi' + . . . , 



E, = MrK ^. = - 2 ^"~''^' ' ^'^ = ~ 2 ^"'-^' ^^ e ^>/o-^'^/,^ 

and M„ , My , M^ , a , and b are determined by the equations 

&hM„ = n;^ — 1 , 

(a + h)M„ — abJf, := — 2 A A , 
Jfo - (a + b) Jl/; + abJI4 = 1 + 7>,2 + 2 A A , 

Jf, - (a + b) Jf, = 2 (AA + A A) , 
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By substituting the values of K„ , A', , and E.^ and multiplying by M^, our 
equation becomes 






,,.,[i_s[»4^]-+.b]i 



,M, 






But 



a + b 



iC, 



AA H 2 ^^•^' ' 






- t, D,M, - - DJJ^D, - />/A - A 



By substituting these, the equation becomes 



a + b 



i/,. 



-a^a+a[ 



1 — 3 



a + bl 






This may easily be transformed into 



z>„^A + a[i + ^.^ + 3A^i ^ + A^ - 1] 



a + bV , 3 



■1 ■ + ^ab 



2-| ' 2"^] 



-aa[i-; 

-|7A[A^ + 2AA-3^^'] 
^a + br3^, 3^.^ S A^A^ -o 
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Which reduces to 



L^ + 2 "~T- ^ + 2 "^T" J ' ' 2 ^' L ^^ - ^ "^ J 

^a + br3^, 3A'An_o 



and thence to 



L ^2 + 2A + l 2A + 1 2A' + lJ ' ~ 

which is perceived to be identical. 
When 

V Q 

is divided by 1 — x the remainder is equivalent to what -= becomes when x 
is put equal to 1. But 

J._ ^ 1{ &—Ie)\x^^) 
y Q Xl-x'-K' ' 

consequently this remainder is 



^ l/(l-a)(l-b) 
^(1) 

the ambiguous sign being so taken as to render the quantity positive. In like 

manner it is shown that the remainder of -7=^ divided by 1 + .« is 

1/^ 



^ •(l + a)(l + b) 
^(-1) 
Then 



^, - + 2 r="5 + ^" + A« + v.^«- + . . . , 

where the upper or lower sign is to be taken according as ^(1) is positive or 
negative. And 

d^, - + 2 1+a; -+ A, + A«+/^2a^ + ---, 
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where the upper or lower sign is to be taken according as H(^ — 1) is positive 
or negative. The expressions for the L and Z', correct to quantities of the 
order of the fourth power of the eccentricities inclusive, are 

2Z„ = [yr^^' -I- z>„ + D, -V a]^; + -^(1) [^; + E.^ , 

2 A =. [if ;^.' - A + a]/^'o + [^^§ + A + a]^; + //(I) E, , 
24 = [a- f^^^ - D, - 2 a] a; + [a- ^f ' - a] ^i + A/^; , 

2z; = [a- ^^' H A + A - a]^, - ^(- 1) [^. - a;] . 

2 A' == [a- ;J|; + A -I a] z^; + [a- ^f " + A + a]^. - //(- 1) E, , 

2Z,' ^ [a- jl^ - A + 2A]^o + [a'^§ + a] ^; + D,E, . 
By substituting the value 



+ b a — b , 

-^- -5j cost's 



and putting 



^, = /, + .,._+i+4[|[i+i]--i.b], 



aV| — ^^, — ^ — y^2 --- 2 — 
^, _ r-a — b _ / . a^ — b^ 



^ab], 



N'. 



_ ,,(a-br 

■> ' '■> o ? 
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where we have, as has been proved above, the relation JV„ = JV^', we get 



y 1 — — 1- - 4- -^^ — cos d) 



-^-- -' ^2 "^^ 



di^y-jOA = :p ^AlJl+M^.^ + iVV + JV,' cos V'' + JV.^ COS 2^'- 4 . . . 



rfc^ ...a + b a — b , 



Integrating, we have 

w — ^6% =c + arc tan [^([-^-5**° ll "^ "^'^ "*" -^i^^"^'' "+ i^2sin2^ + . . . , 

A' -iCt = c'=f arc tan [ J^^t | t^n |1 + iV„Y'' + JV,' sin v'' + i JV,' sin 2^'' + . . . . 

The quadrant in which the arc correspondent to the tangent is to be taken 
is found by dividing the number of the quadrant of </> by 2, if it is even ; or 
by augmenting the number of the quadrant of (/> by unity, if it is odd, and 
then dividing by 2. 

By taking the sine, we have, /9 being any arbitrary angle, 



I ' 1 — X sin {(0 — ^Ci-\- j3) 

= + 1. -T^^ sin ^ cos [W^</' + c + i3 + jy, sin </< + i-lV-i sin ^+ ■■■^ 

+ i,'T^^cos|sin []^„</> + c + ft + iVjSin^/' + i]^,sm2<;' +...], 



l^l +xsm{(o—^Ct + fl) 

= qi ,/r+ asin|'cos [i\^,V + «' + /? + iV^i'sin^'' 4 ^iV^; sin 2^'' + . • •] 

-I- i/i + boos I sin [JVM + e' + /3 + iV,' sin 5/' + JjV; sin 2/' + . . -1 
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or, as they may be written, 



yV^^x sin (co — ^Ct-\- /3) 

=J[|/I=b + v/T=^] sin [(iV;+|)^+6--f-/9+iV,sin9''+ii\^2sin2^''+ . . .] 
+ |[l/r=b ± i/T-a] sin [(iV^, -i)^''+c+;3+iV^, sin ^''+ii\^,sin 2^+ . . .], 



V 1 -h *' sin {(o —\ Ct -[- /9) 

=|[)/TTb + 1/lTa] sin [(iV;,'+|)i.''+6''+;:?+iVi'sini/'+ JiV/sin2i^-h . . .] 
+ 1 [l/T+b ± i/T+li] sin [(i\^;-i)^''+c'+;9-hiV^,'sin^'<+i^2'sin2v''+ . . .]. 

The expression for the auxiliary angle ^'' in terms of the time, which has 
already been obtained, we will denote as follows : 

4' = H. {t + c„) + K, sin 6„ {t + c„) + K, sin 20, {t + e„) + . . . . 

Substituting this for <}' in the preceding formulae, and putting in succession 

;9 = ^Ct, ;9 = 90° + ^Ct, 
we get 

+ P, sin d, {i + c„) + P, sin 2^„ (< + c„) + . . .], 
+ (>, sin d„ {i + 6-,) + Q, sin 2^„ (< + c,) + . . .], 

i/iT^'o° ^y = i[v/T+¥T- i/rq^] ^^^^ [(P„ + i) ^„ {( + .„) + 6'' 

+ /V sin ^„ (< + c„) + P; sin 2d,{t + €,)+...'], 
+ Qi sin d, {t + c„) + $,' sin 2^„ (< + c„) + ...]. 
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Here we have put 

P, --= JV, + (.% + i) A\ , 

P, = \\N, -I- N,K, + 2(i\^, + I) /(-,] , 

Q, = ^1 i\^, + N,K\ -f 2(A^, ~ i) A-l , 

p.: = i-i iv/ + n;k, + 2(iv;, + \) iq , 

It is evident from the equivalent of sin v cos ;- derived from these equa- 
tions that c' = c or c' = c + 180°, according as 

IKh) = l)„ + l\h + n,¥ H- D.p + ... = + vT^l' 
is positive or negative. Hence the latter of the two equations may be written 

1^1 + *ciS ''' = ^ it i/i +~b T i/rq:ii] ^i;^ [(/>„ H- ^) #„ (< 4 '^,) H c 

-f ■ /V sin d„ {( + r„) + /'; sin 2tf„ (< + c,,) ^ . . .] 

± i[ ,/i- -T~b ± i/TT^] f^^ [(P„ - 4) <^o (i -^ Co) + 

+ Q; sin t),{i + c„) + (;/sin2^„(H-c„) + ...], 

where the upper or lower of the newly introduced ambiguous signs is taken 
according as IJ{b) is positive or negative. 



Let us put 



X = {P„ + i) 0„ (t -I- c„) + c, 

i -- {P. - i) ^0 it + f„) + c, 

A = i [l/l^^b + i/l^^, 

A I = i [/l^=~b ± l/i — a], 

A' - + \ [l/l -4- b + l/rTa], 

^\.' -'=- ± h [l/mo ±: |/1 H aj. 
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Then 



+ li AP, + A, (KV + hQ.)] Z (2^ - ^') 
+ [- i A,«. + A (A/V - hP.)\ "j," (2/ - X) 

+ ^{^p;' + ^p.^Z^h-'^t) 

i/^"+^- ci '"' = L A' (1 - yVO + 4 A/<?,'] Z t 

+ [A/(l-i(>,'0 -4A'P,']'jJ/ 

+ [i A'/',' + A,' m:' + 4<;^.')i ci ^^^ - ^'^ 

+ [- kiMQ; + A' (A/-v^ - kPi^ftl^ (2/ - Z) 

+ A,' (A<;>." -*<;>;.) eo" (3;^'- 2z). 

It is evident that e ^ o) and <;' o)' can be expressed in series of the 

cos COS ^ 

same form. 

In applying to Jupiter and Saturn these equations, it is foiind that by 
varying the value of K, 

K'^^ = H 0.0101629, A''^ ^- + 0.0009178, K '^^ = — 0.0050568. 
'iJi dA (IK 

Also 



log [— i/(l — a) (1 — b)J = 9.9574334m, log i (1 + a) (1 + b) = 9.4074864, 

4= -I 0.1672972, A;= I 0.1655301, 

Z, = -I- 0.0028107, LI := — 0.0012760, 

4 = _ 0.0000071, L: = - 0.0000250. 
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Whence 
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Also 



JV„ = + 0.1667632, 
iV, = — 0.0021649, 
K, = — 0.0000021, 

P(, = -f- 0.6837293, 
I\ = — 1434".41, 
P^ = + 5".41, 

e == + 47".17, 



Q-^ 



0".63. 



(i'o 



log A -- 9.5750158, 
log A, r- 0.0101623, 
+ i) d, r... 22".55981, 



iV;; = H 0.1667632, 
iV/ = H 0.0009746, 
iV; = — 0.0000074. 

F; == — 786".82, 
/>' ^ + 2".64, 
Q; = + 694".76, 

<;>; = - 3".50, 

log A' - 9.8634412/*, 
log A,' == 9.7217366, 
(21-^)6,-- 3".50156. 



l/l 



a; ^^° (5 = + 0.3759635 ^^° 

cos COS 

sin 



+ 1.0249824 J^^"^/ 



0.0013085 ^;;; (2,. - x) - 0.0001196 ^^'^ (2,.' - x) 



sin 



+ 0.0000072 ^'j; (3;^ - 2/) H- 0.0000016 ^^^^ (3;.' - 2;^), 



/I — ; — sin -, 
' ' cos 



0.7293089 ®^° y 
cos '• 

0.0013889 ®"^ (2y — /) 
cos ^ '■ '■ ' 

sin 



+ 0.5255160 ^^'^ / 
' cos '■ 



0.0008842 '"" (2/ — y) 
cos ^ '• '-' 

, sin 



- 0.0000058 '^^^^ (3y - 2/) + 0.0000052 ^^^^^ (3/ - 2y). 



The value of c is found to be 

c = 340° 8' 50".26. 

Hence the expressions for the two arguments are 

y = 308° 13' 15".13 -| 22".55981i!, 
y' =. 31° 4' 5".98 ^- 3".50156<. 
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The following expressions for e and e' were obtained : 

[8.6282138] ,/(l — [6.5410419] cos ^''), 
[8.8231642] ^/(l + [6.9312571] cos ^'0, 



e 



1 1 — X 



V T — X 

e' 
v'l -\- X 

= [8.6282135] [1 — [6.24001] cos(x — /) + [3.7900] cos2(x — /)}, 



e 



= [8.8231648] {1 + [6.63023] cos {^ — /) — [4.1982] cos 2 (j^ — /)}. 



1 1+3- 

By means of these we can pass to the expressions for the following functions : 

e ^^^ M r--= + 0.01596822 ®^° y + 0.04354278 ^^^ / 

cos COS '• cos '- 

— 0.00005696 ^*" (2y — /) — 0.00000886 ^^° (2/ — y) 

COS ' ' COS " 

+ 0.00000031 ®^° (3y — 2y') + 0.00000009 ^^^ (3/ — 2y), 

COS ' ' COS 

, sin -, rv njornnriA sin , n no A nc A t\n ^Va 



e CO 



— 0.04852990 "'" y + 0.03496407 '''" / 



COS cos '• COS 

^- 0.00008205 ^^^ (2y — /) — 0.00005134 ^^^ (2/ — y) 

COS ^ '• '• ' COS ^ '' '•' 

-^ 0.00000033 ®^° (3y — 2y') + 0.00000031 ^'^^ (3/ — 2y). 

COS ' COS " 

It will be observed that these expressions are as convergent as could be 
wished. The form of these integrals being discovered, another and more 
direct method of arriving at them is suggested. The coefficients being 
assumed as indeterminate as well as the rates of movement of the two 
arguments together with the constants which complete the values of the 
latter, the expressions could be substituted in the differential eqiiations, and 
thus would arise twelve equations of condition, which along with the value 
of the four variables at the origin of time would determine the sixteen 
unknowns involved. But on trial it seems this way of proceeding would 
necessitate as long computations as the method we have followed. 

In conclusion, it may be observed that, if terms arising from the squares 
and higher powers of the masses were taken into consideration, the form of 
this investigation would not thereby be changed ; the only effect produced 
would be that the values of the various constants involved would receive 
slight modifications. 



